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ABSTRACT 

Methods  of  computation  and  error  estimates  are  presented  for  numerical 
appi oxi mation  of  selected  meteorological  parameters  used  in  the 
solution  of  cloud  physics  problems  These  parameters  are  latent 
heat,  saturation  vapor  pressure,  dew  point  and  wet  bulb  temperatures, 
specific  heat  of  liquid  water,  temperature  of  fhe  lifting  condensation 
level,  change  in  mass  of  #realy  falling  water  drops,  and  saturation- 
adiabalic  temperatures.  The  FORTRAN  programming  language  listings 
of  the  approximations  are  provided  in  the  Appendix. 


\ 


5 


i,i<  '**•' 


CONTENTS 


1 


INTRODUCTION  . 

METEOROLOGICAL  PARAMETERS  . 

Latent  Heats  . 

Saturation  Vapor  Pressures  . 

Dew  Point  Temperature . 

Wet  Bulb  Temperature  . 

Specific  Heat  of  Liquid  Water  . 

METEOROLOGICAL  PROCESS  PARAMETERS  . 

Temperature  of  the  LCL . •  • 

Change  in  Mass  of  Freely  Falling  Water  Drops  .  . 

Saturation-Adiabatic  Temperatures  . 

CONCLUDING  REMARKS  . 

LITERATURE  CITED  . 

APPENDIX  —  COMPUTER  LANGUAGE  LISTINGS  OF  FUNCTIONS 

AND  SUBROUTINES  . 


Page 


2 

2 

4 

7 

10 

12 

12 

12 

14 

17 


24 


25 


27 


v 


LIST  OF  TABLES 


TEMPERATURE  UNIT  CONVERSIONS 


Teal  g  UNIT  CONVERSIONS 


PERCENTAGE  ERRORS  IN  THE  NUMERICAL  APPROXIMATIONS 
FOR  THE  CA-CULATION  OF  LATENT  HEATS  . 


PRESSURE  UNIT  CONVERSIONS 


PERCENTAGE  ERRORS  IN  THE  NUMERICAL  APPROXIMATIONS 
FOR  THE  CALCULATION  OF  SATURATION  VAPOR  PRESSURES  AT 
SELECTED  TEMPERATURES  . 


PERCENTAGE  ERROR  OF  THE  DEW  POINT  aND  WET  BULB 
TEMPER/', TURE  ITERATION  TECHNIQUES . 


PERCENTAGE  ERRORS  IN  THE  NUMERICAL  APPROXIMATIONS  FOR 
THE  CALCULATION  OF  THE  SPECIFIC  HEAT  OF  LIQUID  WATER  . 

ERROR  (C)  OF  THE  LIFTING  CONDENSATION  LEVEL  TEMPERATURE 
FUNCTION  . 

PERCENTAGE  ERROR  IN  THE  NUMERICAL  APPROXIMATIONS  FOR 
THE  CALCULATION  OF  dM/dt  AT  SELECTED  OONDIT’ONS.  .  . 


COMPARISON  OF  THE  RESULTS  OF  THE  NUMERICAL  METHOD  FOR 
SOLVING,A  DIFFERENTIAL  EQUATION  APPLIED  TO  dx/dy  = 

+1 . 


XI.  COMPARISONS  OF  SATURATION-ADIABATIC  TEMPERATURES 


22 


Symbol  Subscri  pt 
c 

w 

e 

3 

a 

s 

s. 


log 

L 


M 

P 

R 

rh 

t 

T 


a 

t 


a 

d 

LCL 

m 

s 

w 

x 


s 

w 

z 


LIST  OF  SYMBOLS 
Item 


Sped  f  ic  Heat 

air  at  constart  pressure 
I iquib  water 


Units  Used 
ITca I  g  *  K 


Water  Drop  Diameter  m 

Partial  Pressure  of  Water  Vapor  mb 

actual 
saturation 

saturation  over  sol  id  water 
saturation  over  liquid  water 
saturation  at  the  wet  bulb  temperature 

Common  Logarithm  (base  10) 

Latent  Heat  of  water  I  Tea  I  g 

fus ion 
subl imation 
vaporization 


Mass  of  Water  Drop 
Pressure 

atmosphere  (with  no  water) 
atmosphere  total 

Dry  Air  Gas  Constant 

Re  I  at ive  Humidity 

T  i  me 

Temperature 
dry  bulb 
dew  po i nt 

lifting  condensation  level 
psue do-adiabatic 
saturat ion-ad i abat ic 
wet  bul b 

approximation  of  saturation- 
adiabatic 

Mixing  Ratio 

water  vapor  at  saturation 
I  iqu id  water 
initial  val ue  of  X 

s 

Ratio  of  molecular  weights  of 
water  vapor  and  dry  air 


-I 


9 

mb 

I  Tea  I  g  * K 
Pe  rcent 
Sec 
K 


9  9 


VI  i 


INTRODUCTION 

The  application  of  techniques  for  numerical  simulation  to 
various  meteorological  processes  requires  that  many  rreteoro  logica  I 
parameter  descriptions  be  subject  to  individual  approximation. 

Thus,  the  accuracy  and  adequacy  ot  the  descriptions  are 
dependent  on  the  requirements  of  the  individual  "simulator." 

In  many  cases  the  importance  of  the  meteorological  parameter 
description  is  subordinate  to  some  broader  scale  problem;  as 
such,  if  the  description  is  deemed  adequate,  it  is  normally 
used  without  explanation  or  comment.  As  a  consequence 
there  are  many  methods  available  for  the  description  of 
meteorological  parameters,  but  there  is  a  paucity  of  information 
on  how  well  they  describe  the  parameters.  If  the  accuracy 
must  be  known,  then  the  investigator  must  either  construct 
his  own  approximation  or  ascertain  the  accuracy  of  some 
available  roproxi mation,  either  of  which  can  detract  from 
the  solutio,"  of  a  larger  problem. 

This  report  contaii.s  descriptions  of  numerical  approximating 
techniques  for  meteorological  parameters  required  by  the  author 
and  his  associates  in  conjunction  with  numerical  modeling  efforts 
concerning  clouds  and  fogs  [l],  [23,  [33,  [43.  This  presentation 
is  made  so  that  there  m!ght  be  a  ready  reference  to  a  set 
of  approximating  methods  and  the  accuracies  associated  with 
the  calculated  parameter  values. 

The  numerical  approximations  are  divided  into  two  groups: 

(I)  descriptions  of  meteorological  parameters,  per  se;  and  (2) 
descriptions  of  rreteoro  logical  parameters  associated  with  a 
meteorological  process.  The  firsT  group  covers  the  latent 
heats  of  water,  the  saturation  vapor  pressures  of  water,  the 
dew  point  and  wet  bulb  temperatures,  and  the  specific  heat  of 
I iquid  water.  The  second  group  covers  the  temperature  of  the 
lifting  condensation  level,  the  change  of  mass  of  freely 
falling  water  drops,  and  the  temperatures  resulting  from 
saturated  adiabatic  processes.  The  method  of  approximation 
and  the  source  of  comparative  values,  as  well  as  the  error 
as  determined  relative  to  the  comparative  values,  is  provided  for 
each  parameter.  The  units  for  the  comparisons  are  those  of  the 
comparative  base,  where  applicable.  In  those  cases  where  comparative 
values  are  not  available,  the  basis  of  the  accuracy  given  is  discussed. 
When  percentage  errors  are  given  the  value  presented  is  calculated 
from 


Percentage  Error  =  100 


Calculated  val ue -Reference  value 
Reference  v3iue 


(I) 


The  FORTRAN  programming  language  listing  of  the  approximations 
is  provided  in  the  Appendix  in  either  FUNCTION  STATEMENT  o'" 
SUBROUTINE  form  as  aporopriate. 


METEOROLOGICAL  PARAMETERS 

The  rreteoro logical  parameters  described  in  this  section  are 
associated  witn  the  state  of  the  atmosphere  and  are,  therefore, 
essentially  definitions  of  the  value  of  the  parameter  for 
the  specified  state.  The  parameters  covered  are  those  associated 
with  the  following:  the  heat  requ i rements  for  the  change  of  phase 
of  water;  the  partial  pressure  of  water  vapor  over  liquid  and 
sol  id  plane  water  surfaces;  the  Temperature  at  saturation  of  air 
ccoled  at  constant  pressure  with  no  change  in  water  vapor  content; 
the  temperature  of  air  cooled  to  saturation  at  constant  pressure; 
and  the  specific  heat  of  liquid  water. 


Latent  Heats 


The  latent  heats  are  associated  with  the  vapor-liquid,  vapor-solid, 
and  liquid-solid  phase  changes  of  water.  The  latent  heats  can  be 
expressed  as  functions  of  temperature,  T,  with  the  techniques 
given  here  using  degrees  Kelvin  (conversion  to  other  femperature 
units  is  provided  in  Table  I)  and  the  latent  heat  expressed 
in  ITcal  g  (conversion  toother  units  is  provided  in  Table  II). 

The  latent  heat  associated  with  the  liquid-vapor/vapor-liquid 
phase  changes  is  the  latent  heat  of  vaporization,  L  .  It 
can  be  approximated  from 

L  =  754.817  -0.575  T  (2) 

v 


The  latent  heat  associated  with  the  solid-vapor/vapor-solid  phase 
changes  is  the  latent  heat  of  sublimation,  L  .  It  can  be 
approximated  from 
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TABLE  I1 

TE^ERATURE  UNIT  CONVERSIONS 


Unit 

Symbol 

Unit 

Symbo  1 

Conversion 

Kelvin 

K 

Centigrade 

C 

K  -  273.16 

Fahrenheit 

F 

1.8  (K  -  273. 16)  + 

32 

Rankine 

1.8  (K  -  273. 16)  + 

491 .69 

TABLE  II2 

ITcal  g'1  UNIT  CONVERSIONS 

1.0 

Internationa  1 

Steam  Tables 

Calorie  (ITcal)  =  1.00032  cal^ 

LG  ITcal  g  *  =  4.18684  x  10^  cm2  sec  2 
4. 1 8684  x  I03  m2  sec  2 
4.5057  x  I04  ft2  sec"2 


The  conversions  given  here  are  based  on  the  values  provided  by 
List  [5,  p.  173. 

^The  conversions  given  here  are  based  on  the  values  provided  by 
List  [5,  p.  14]. 
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(3) 


l_s  =  98.777  +  [335566.953  -  |  T  -238|  2D  ^ 

The  latent  heat  associated  with  the  solid-liquid/liquid-solid 
phase  changes  is  the  latent  heat  of  fusion.  L^.  It  can  be 
approximated  from 

Lt  =  (1.7757  -0. 002324 5T )T  -  0.012  |T-243|  -231.494  (4) 

i 

The  accuracy  of  Lv,  Lg,  and  is  given  in  Table  III. 


Saturation  Vapor  Pressures 


In  neteorological  parlance,  saturation  vapor  pressure  refers  to 
the  partial  pressure  of  water  vapor  for  saturated  conditions  over 
a  plane  surface  of  either  liquid  or  solid  water.  The  saturation 
vapor  pressure,  expressed  in  millibars  (conversion  to  other 
pressure  units  is  provided  in  Table  IV),  is  a  function  of  temperature 
only. 

The  saturation  vapor  pressure  associated  with  a  plane  liquid 
water  surface,  e  ,  is  approximated  by 
w 

. „(a+b+c+d+e)  /c. 

e  =  I  u  (o ; 

s 

V 


where  a 
b 
c 
d 
e 


23.83224 
-5.02808  log,T 

i  r\ 


8- 1  328 10  ' 
-1.3816  x  10 
-2949.076/T 


x  10 


(3.49 f-J?  -  I  302.8844/T) 


x  10 


(11.344  -  0.0303998  T) 


and  the  saturation  vapor  pressure  associated  with  a  plane  solid 
water  surface,  eg  ,  is  approximated  b; 


Equations  5  and  6  are  restaterrenl s  of  those  given  by  List 
[5,  p.  350]. 


I 


TABLE  II r 

PERCENTAGE  bRPORS  IN  THE  NUMERICAL  APPROXIMATIONS  FOR  CALCULATION 

OF  LATENT  HEATS 


Percenlage  Errors  in  Latent  Heats 


Temperature 

(C) 


Vaporization  Sublimation  Fusion 


60 

0.009 

55 

0.004 

50 

0.000 

45 

-0.004 

40 

0.009 

35 

0.004 

30 

0.017 

25 

0.030 

20 

0.T  43 

15 

0.C/.8 

10 

0.051 

5 

0.0;.  ; 

0 

0.074 

0.000 

0.063 

-10 

0.083 

0.003 

0.108 

-20 

0.057 

-0.005 

-0.080 

-30 

0.000 

0.005 

-0.246 

-40 

-0.153 

0.005 

-0.462 

-50 

-0.445 

-0.005 

0.350 

-60 

0.003 

-70 

0.000 

-80 

0.001 

-90 

0.006 

-100 

0.000 

The  comparative  base  for  construction  of  this  Table  was  that 
provided  by  List  \J>,  p.  3433. 


TABLE  IV5 

PRESSURE  CONVERSIONS 


3  -I  -2 

1 .0  mi  i I ibars  (mb)  =10  g  cm  sec 

1.4504  x  I0"2  lb  in”2 


2.08854  lb  ft 


-2 


TABLE  V6 

PERCENTAGE  ERRORS  IN  THE  NUMERICAL  APPROXIMATIONS  FOR  CALCULATION 
OF  SATURATION  VAPOR  PRESSURE  AT  SELECTED  TEMPERATURES 


Ftercen+age  Errors 


Temperature 


(C) 

e 

s 

w 

e 

s  j 

-;oo 

0.003 

-80 

-0.015 

-60 

0.031 

-40 

0.017 

0.008 

-20 

-0.003 

-0.036 

0 

-0.004 

-0.008 

20 

-0.004 

40 

-0.003 

60 

-0.004 

80 

-0.005 

100 

-0.004 

The  conversions  given  here  are  based  on  the  values  provided  by 
l ist  [5,  p.  Ijj. 

^Ihe  comparative  base  for  the  construction  of  this  table  was  that 
provided  by  List  [5,  pp.  351-353  and  360-3613- 
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>  _  |Q<a+b+c+d) 

s. 

! 


(6) 


where  a 
b 
c 

and  d 


2.0702 
-2484. 99/T 
3.56654  log  T 
-0.0032098  T 


The  percentage  accuracies  of  the  functions  of  e  and  e 
are  given  in  Table  V^.  Sw  s 


Dew  Point  Temperature 


The  dew  point  temperature,  T  ,  is  that  which  unsaturated  air,  cooled 
at  constant  pressure  with  no  change  in  the  water  vapor  mixing  ratio, 
would  have  at  saturation.  The  process  for  approximating  this  temperature 
is  an  iterative  one,  with  the  basis  of  closure  being  the  difference 
between  the  actuai  vapor  pressure  determined  from  the  dry  bulb 
temperature  and  relative  humidity  and  the  calculated  vapor  pressure 
of  the  dew  point  temperature.  Assuming  that  the  dry  bulb  temperature 
and  the  relative  humidity  are  known,  the  calculations  are  accomplished 
as  foi lows: 


I.  Select  the  accuracy  desired  in  the  vapor  pressure 

differences.  The  resu It i ng  number  of  iterations  required  for 
various  accuracies  are  provicad  in  Table  VI. 


Calculate  the  saturation  vapor  pressure  of  the  dry 
bulb  temperature,  e$,  using  Eq.  5. 


3. 


Calculate  the  desired  vapor  pressure  at  the  dew  point 
temperature,  e  ,  from  the  relationship 
e  rh  3 


e 

a 


s 

100 


(7) 


4.  Make  an  initial  estimate  of  the  dew  point  tenperature. 
The  initial  estimate  in  this  report  is  found  from  a 
restatement  of  an  approximation  given  by  -lew^n  and 
Longley  C6,  p.  2263: 


Applying  the  function  statements  stated  here  to  all  the  values 
providea  by  List  [5,  pp.  351-353;  360-3613  the  maximum 
percentage  error  found  was  0.04  percent. 


2540 


(8) 


8.573  -  !og(O.I6367ea> 


5.  Establish  an  initial  increment  for  tenperature  iteration. 
The  initial  increment  of  tenperature,  AT,  used  in  this 
report  is  determined  from 


AT  =  7.7  -0.026T 


(9) 


with  the  added  condition  that  |at|>o. I . 

6.  Calculate  the  saturation  vapor  pressure  corresponding 

to  the  estimated  T^,  using  Eq.  5  of  Table  V. 

7.  Find  the  error  in  this  vapor  pressure  from 


error  = 


100  Te  -  e  (T  )} 
a  s  d 


(10) 


If  this  error  is  smaller  than  the  value  specified  in 
Step  I,  STOP:  the  T^  currently  being  investigated  is 
acceptable  as  an  estimate  of  the  dew  point  tenperature. 

8.  If  the  sign  of  the  error  has  changed  from  tie  previous  one 
and  the  error  is  larger  than  the  value  specified  in  Step  I, 
change  At  to  *sAT. 


9.  Adjust  the  estimated  value,  T^,  according  to 


=  Td  + 


e  -  e  CT  ,) 
a  s  d 

je  -  e  (T  ) | 
a  s  d 


AT 


(II) 


and  then  RETURN  TO  STEP  6. 

The  accuracy  of  the  dew  point  tenperature  is  dependent  on  the 
accuracy  of  Eq.  5  and  the  temperature  of  the  dew  point.  With  a 
closure  accuracy  requirement  of  0.001  percent,  the  maximum 
error  to  be  expected  would  range  from  0.005  C  at  -50  C  to 
0.002  C  at  +100  C. 


The  percentage  accuracies,  exclusive  of  consideration  of  the 
vapor  pressure  function  contributions,  are  provided  in  Table 
Vi . 
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Wet  Bulb  Temperarure 


The  wet  bulb  temperature,  T  ,  is  that  which  ynsaturated  air,  cooled 
by  the  adiabatic  evaporation  of  water  at  constant  pressure  with 
no  external  source  of  heat,  would  have  at  saturation.  Like  the 
dew  point  temperature  of  the  preceding  section,  the  process  for 
approximating  this  temperature  is  ar.  iterative  one,  with  the  basis 
of  closure  being  the  difference  in  vapor  pressures.  Assuming  that 
the  dry  bulb  temperature  and  relative  humidity  are  known,  the 
calculations  are  accomplished  as  follows: 

1.  Select  the  accuracy  desired  in  the  vapor  pressure 
differences.  The  resulting  number  of  iterations  required 
for  various  accuracies  are  provided  in  Table  VI. 

2.  Calculate  the  saturation  vapor  pressure  of  the  dry  bulb 
temperature,  e ,  using  Eq.  5. 

3.  Calculate  the  vapor  pressure  at  the  dew  point  temperature, 
e^,  from  the  relationship  stated  in  Eq.  7. 

4.  Hake  an  initial  estimate  of  the  wet  bulb  temperature. 

The  initial  estimate  in  this  report  is  foimd  from  a 
consideration  of  the  dew  point  tenperature,  T  ,  as 
estimated  in  Eq.  8,  in  +he  form 

e  T  +  e  T  . 

T  =  -  ~  -  ■ - —  (12) 

w  e  +  e 

a  s 


5.  Make  an  initial  estimate  of  the  saturation  vapor  pressure 
of  the  dew  point  temperature,  e  ,  using  tne  expression 
from  List  [ '5 ,  p.  3661  in  the  form 


e  =  e  +  0.000(3  P(l  +  0.001 1  T  )(T  -T  )  (13) 

w  a  w  a  w 


6.  Establish  an  initial  increment  for  temperature  iteration. 
The  initial  increment  of  temperature,  AT,  used  in  this 
report  is  determined  from 


AT  »  [e  (T  )  -  e  ]  (3.3  -  0.0089T  ) 

s  w  w  a 


(14) 


where  e„(T  )  ss  found  from  Eq.  5,  with  the  added 
condition  ¥hat  j dT ! >0-2- 


7.  Calculate  the  saturation  vapor  pressure  corresponding  to 
the  estimated  T^,  e  (T^),  using  Eq.  5. 

8.  Calculate  the  vapor  pressure,  e^,  using  Eq.  13. 

9.  Find  the  error  between  the  two  vapor  pressures  (from 
Steps  7  and  8)  from 


error  = 


[e  (T  )  -  e  ]  100 

s  w  w 


(15) 


If  th;s  error  is  smaller  than  the  value  spec’ tied  in 
Step  I,  STOP:  the  T  currently  being  investigated  <s 
acceptable  as  an  estimate  of  the  wet  bulb  temperature. 

10.  If  the  sign  of  the  error  has  changed  from  the  previous  one  and 
the  error  is  larger  than  the  value  specified  in  Step  I,  change 
AT  to  *sAT. 

11.  Adjust  the  esfinerea  value,  T  ,  according  1o 


T  =  T  +  AT 

w  w 


(16) 


where  AT  is  positive  i f  e  (T  )  <  e  and  negative  if 

e  (T  )  >  e  .  RETURN  TO  §TEP  7.  w 

S  w  i 

As  for  the  dew  point  temperature,  the  accuracy  of  the  wet 
bulb  tempera+,Jr"e  is  mainly  dependent  on  the  accuracy  of  Eq.  5 
and  the  tempera; ur°  of  the  wet  bulb. 
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Specific  Heat  of  Liquid  Water 


The  amount  of  heat  required  to  raise  one  gram  of  liquid  water  one 
degree  K  is  the  specific  heat  of  water,  c  .  The  function  for  estimating 
the  value  of  c  results  from  a  curve  fitting  cf  the  data  given 
by  List  C5,  p.*343l.  The  function  is 


c 

w 


0.9979  + 


0.3021 

2 (8.5—0. f  | T-308 | ) 


(17) 


The  accuracy  of  the  application  of  this  function  with  the 
values  from  List  £5,  p.  343]  is  summarized  in  Table  VII. 


METEOROLOGICAL  PROCESS  PARAMETERS 

The  meteorological  parameters  described  in  this  section  are  associated 
with  meteorological  processes.  The  parameters  covered  are:  (I) 
the  temperature  air  would  have  after  it  has  been  lifted  (dry  adiabatical ly 
with  no  change  in  water  vapor  mixing  ratio)  from  an  initial  level 
until  saturation  is  reached,  which  corresponds  to  the  temperature 
of  the  lifting  condensation  level  (LCL);  (2)  the  time  rate  of  change 
cf  the  mass  of  freely  falling  water  drops;  and  (3)  the  temperature 
resulting  from  moving  air  from  one  pressure  level  to  another  in 
a  maintained  saturated  state,  which  is  the  psuedo-adiabatic  temperature. 


Temperature  of  the  LCL 

The  temperature  of  the  LCL,  T  ,  is  computed  from  the  dry 
bulb  temperature,  T,  and  the  dew  point  temperature,  Trf,  of 
the  initial  level  in  a  restatement  of  the  formulation  provided 
by  Inman  £7],  namely 


T,  -  T ,  +  (0.098  -0.001571  T,  +  0.000436  T  )  (T  -T  )  (18) 

LCL  d  d  a  a  d 
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TABLE  VII 


PERCENTAGE  ERROR  IN  THE  NUMERICAL  APPROXIMATIONS  TOR  THE 
CALCULATION  OF  THE  SPECIFIC  HEAT  OF  LIQUID  WATER 


Temperature 


Percentage  Error  of 


(C) 

c 

w 

-50 

0.000 

-40 

0.785 

-30 

-0.609 

-20 

-0.417 

-1C 

-0.316 

C 

-0.006 

5 

0.087 

10 

0.132 

15 

0.144 

20 

0.146 

25 

0.127 

30 

0.108 

35 

0.084 

40 

0.108 

45 

0.137 

50 

0.176 

55 

0.234 

60 

0,322 

13 
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The  accuracy  of  the  approximation  is  given  in  Table  V i II ,  where 
the  comparative  base  is  that  provided  by  Nordquist  £8]. 

These  values  were  computed  using  an  iterative  scheme,  which 
had  a  closure  of  0.001  percent  in  determination  of  the  LCL  water 
vapor  mixing  ratio. 


Change  of  Mass  of  Freely  Falling  Water  Props 

The  construction  of  numerical  models  which  contain  consideration 
of  liquid  water  require  that  some  method  be  established  for 
considering  the  mechanism  of  evaporation.  The  parameterization 
used  here  is  based  on  the  method  proposed  by  Kessler  [9,  pp. 
29-303  using  the  data  provided  by  Kinzer  and  Gunn  DOj. 

The  basis  of  the  approximation  is  that  the  change  in  mass  due  to 
evaporation  can  be  expressed  in  the  form 


~  =  A  Ce  -e  )  dB  (19) 

dT  s  a 

The  two  approximations  provided  here  are  based  on  "best" 
fits  of  the  Kinzer  and  Gunn  data  to  Eq.  (19),  one  as  the 
equation  stands  and  the  other  including  a  tenperature 
dependence.  The  first  approximation  is 


=  0.0152  (e  -e  )  d  1  *615  (20) 

dT/j  s  a 


while  the  second  is 


f)2  =[l.  -  0.025(T,-290>]  $),  (21. 

The  accuracy  of  these  approximations  for  selected  conditions  is 
given  in  Table  IX,  with  the  Kinzer  and  Gunn  data  as  the 
comparative  base. 
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TABLE  VIII 

ERROR  (C)  OF  THE  LIFTING  CONDENSATION  LEVEL 
TEMPERATURE  FUNCTION 


TABLE  IX 


10 


PERCENTAGE  ERROR  OF  THE  dM/dt  APPROXIMATIONS  FOR  FREEI  Y 
FALLING  WATER  DROPS  FOR  SELECTED  CONDITIONS 


_ Temperature _ 

Drop  Relative 

Diameter  0  C  40  C  Humidity 

(cm)  (%) 


(dM/dt) 

(dM/dt ) 

(ar7dt ) 

(dM/dt) 

0.02 

50 

29 

-43 

40 

10 

45 

29 

-84 

38 

50 

44 

28 

-90 

36 

90 

0.08 

33 

4 

-96 

17 

10 

33 

5 

-102 

15 

50 

32 

3 

-109 

12 

90 

0.20 

32 

3 

-106 

13 

10 

32 

4 

-112 

1  1 

50 

31 

2 

-!  19 

8 

90 

0.30 

43 

18 

-84 

22 

10 

43 

19 

-89 

20 

50 

42 

17 

-96 

17 

90 

0.36 

50 

29 

-64 

31 

10 

50 

29 

-63 

29 

50 

49 

28 

-'4 

27 

90 

While  these  percentage  errors  are  quite  large,  they  relate  to  rather 
small  numbers,  and  are  adequate  for  those  problems  where  the  cloud 
physics  are  treated  in  a  gross  sense. 


Saturation-Adiabatic  Temperatures 


Air  roved  adiabatical I y  from  one  pressure  level  to  another  in  a 

saturated  state  can  take  on  one  of  two  different  values, 

depending  on  the  disposition  of  the  water  content  of  the  air. 

If  the  total  water  content  remains  constant,  the  process  is 

reversible  and  the  resulting  temperature  is  the  saturation-adiabatic 

temperature,  T  .  If  the  total  water  content  decreases  by  the  amount 

of  liquid  water  present  s'  that  the  process  is  irreversible, 

the  resulting  temperature  is  the  psuedo-adiabatic  Temperature, 

T  . 
m 


The  technique  for  determining  these  temperatures  requires  use 
of  the  statemenf  of  the  First  Law  of  Thermodynamics,  based  on 
the  definitions  provided  by  Beers  [II,  p.  359],  in  the  form 


c  dT 
P 


RT 


(pt-v 


L  X 

d(P  -e  )  +  T  d(  V  -> 

t  s  T 


+  c  (X  +  X  )dT  =  0 
w  s  w 


(22) 


where 


€  e 

(FviT) 


(23) 


and  the  Clausius-CI apeyron  Equation  is 


de 


KL  e 
v  s 


RT2 


dT 


(24) 


Differentiation  of  Equations  (2)  and  (23),  and  substitution 
along  with  Equation  (24)  into  Equation  (22),  yields  after 
rearrangement 
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dT 

<P 


(P.-e  ) 
t  s 


+ 


RT  +  L  X 
v  s 


c  (X  +  X  >-0.575X  + 

W  S  w  S 


TT 
X  L  P 


s  v 


RT2e 


(25) 


To  establish  a  value  for  the  saturation-adiabatic  tenperature, 
it  is  assumed  That  no  hydroneteors  are  within  the  air  at  the 

start  of  the  process,  so  that  (X  +  X)  =  X  (T  ,  P  )  =  X  . 
Thus  s  w  s  o  o  z 


RT  +  L  X 

_ S _ V  s 

(P  -e  Me  +  c  X  -0.575X  + 
t  s  rp  w  z  s 


x2l2p. 

s  v  t 
RT2e 

s 


(26) 


For  the  psuedc-adiabatic  tenperature,  it 
the  liquid  water  is  removed  from  the  air. 
Thus 


s  assumed  that  a  1 1 

so  that  X  =0. 

w 


RT  +  L  X 
m  vs 


Wp7 

S  V  t 

RT2e 


(P.-  e  i  +  c  X  -0. 
t  ->p 


..e  -575X  + 

w  s  s 


(27) 


For  many  meteorologic? I  problems  it  is  sufficient  to  assume  that 
the  partial  pressure  exerted  by  thr  air  is  much  larger  than  that  exerted 
by  the  water  vapor  included  in  the  air,  that  the  relative  mass  of  the 
water  vapor  is  smal I  compared  to  the  mass  of  the  air,  and  that  the 
latent  heat  of  vaporization  is  constant.  With  these  assumptions,  Eq. 
(22)  can  be  expressed  as 


c  dT  -  -  dP,  t  L  dX  -=0  (28) 

p  P+  t  vs 
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and  £q.  (23)  can  be  expressed  as 


(29) 


He:  e  T  is  used  to  designate  this  aeteoro logical  approximation 
to  thexpsuedo-adiabatic  tei^perature 


dT 


dP, 


RT  +  L  x 
x  e  s 


X2L2P 
s  v  t 

RT2« 


(30) 


where  for  each  Equation  (26),  (27),  and  (30),  e  and  L  are 

assumed  to  be  functions  of  T  ,  T  ,  or  T  ,  as  appropriate. 

s  m  x  vy  y 

The  solution  of  these  equations  is  sought  by  using  a  second-order 
predictor-corrector  method  Cl2,  pp.  385-386],  where  fhe  initial 
step  is  accomplished  using 


Thereafter,  a  first  approximation  "predictor"  is  made  using 


AP 


(32) 


and  the  "corrector''  is 


1  (A  systems  sub-program  for  the  Univac  1108  computer  is  available 
for  accomplishing  this  integration  in  a  more  general  manner 
Cl 8,  sec.  I  I ,  p. I ]. 
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€55? 


The  "corrector’*  is  repeated  ujtil 


specified  closure  C34) 


To  provide  some  indication  of  the  validity  o*  this  integration 
technique,  the  egression 


dx 

dy 


0.5 


(35) 


was  integrated  froqa  y=l  to  y=30,  with  a  ay  of  0.01  and  a  closure 
requirement  of  10  .  The  results  are  given  in  Table  X. 

While  an  accuracy  analysis  of  the  saturation-  and  psuedo-adi abatic 
temperatures  is  not  possible  due  to  a  lack  of  a  suitable 
comparative  base,  it  is  possible  to  compare  the  results  with  those 
obtained  by  another  approximation  method  C5,  p.  3183.  The 
comparisons  are  given  in  Table  XI  and  include  the  values  provided 
by  List  C5,  pp.  320-3213  for  psuedo-adi abatic  temperatures 
with  initial  conditions  of  1096.7  mb  and  I4C,  program 
requirements  of  Ap  =  2  mb,  and  a  closure  requirement  of  10  . 
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TABLE  X 


€£»«WI!S£»i  OF  TiH£  RESULTS  OF  THE  NUGRICAIL  MEBC&  FCR  SOLVING  A 
OlFFSROtfTlAL  EQUATION  APPLIED  TO  «tx/«*y  =  %  +  ! 


Value  of  y 

Initial 

of 

Value 

V 

Value  of 

Integration 

X 

% 

x  =  y  +  y 

Difference  of 

Values  of  X 

Percentage 

2 

1 

3.41421 

3.41421 

-0.00004 

4 

1 

6.00000 

6.00000 

-0.00006 

6 

1 

8.44949 

8.44949 

-0.00006 

8 

1 

10.82843 

10.82843 

-0.00006 

16 

1 

20.00000 

20.00001 

-0.00005 

25 

1 

30.00000 

30.00001 

-0.00004 

30 

1 

35.47723 

35.47723 

-0.00002 

1/&LE  XI 


€0#f%£8sms  Of  SSTllRftT !C8S-A3$BASA7 SC  TElf’ERftl&iBES  CO 


PRESSURE  PS£€G-ffllfimTIC  smmTE®-ffll*3AJlC 


LIST 

CaiOTtATED 

C*mHLATED 

J€T  /FFWSX 

1096.7 

14.0 

14.0 

14.0 

14.0 

5046.7 

12.0 

12.0 

12.0 

12.0 

1000.0 

10.0 

10. 1 

10.0 

10.0 

957.0 

8.0 

8.1 

S.l 

8.1 

915.0 

6.0 

6.1 

6.1 

6.1 

876.0 

4.0 

4.2 

4.1 

4.1 

840.0 

2.0 

2.2 

2.2 

2.1 

805.0 

.0 

.2 

.1 

.1 

773.0 

-2.0 

-1.7 

-1.8 

-1.9 

744.0 

-4.0 

-3.6 

-3.7 

-3.8 

715.0 

-6-0 

-5.6 

-5.7 

-5.8 

688.0 

-8.0 

-7.5 

-7.7 

-7.8 

663.0 

-10.0 

-9.4 

-9.6 

-9.7 

638.0 

-12.0 

-11.4 

-11.7 

-11.8 

616.0 

-14.0 

-13.3 

-13.6 

-13.7 

595.0 

-16.0 

-15.2 

-15.5 

-15.6 

573.0 

-18.0 

-17.2 

-17.6 

-17.8 

554.0 

-20.0 

-19. 1 

-19.6 

-19.7 

536.0 

-22.0 

-20.9 

-21.5 

-21.6 

518.0 

-24.0 

-22.9 

-23.5 

-23.6 

502.0 

-26.0 

-24.7 

-25.3 

-25.5 

485.0 

-28.0 

-26.7 

-27.4 

-27.6 

470.0 

-30.0 

-28.5 

-29.3 

-29.5 

454.0 

-32.0 

-30.6 

-3!  .4 

-31.6 

440.0 

-34.0 

-32.4 

-33.4 

-33.5 

426.0 

-36.0 

-34.4 

-35.4 

-35.5 

413.0 

-38.0 

-36.2 

-37.3 

-37.5 

400.0 

-40.0 

-38.1 

-39.3 

-39.4 

387.5 

-42.0 

-40.0 

-41 .3 

-4|.4 
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r 

h. 


las l£  Xt  iosts'i 


memssaus  o?  sftTni!mTfi9iiH®B«afti!:c  faifJss*B(iiRes  «o 


P5€SSiia£  P5BI0BHttMA3fcT6e  SJfflW^aiHMMiaKTiC 


LIST 

CaiOJSATE® 

MET  #P!P80!3t 

375.3 

-44.0 

-42-0 

-43.3 

-43.4 

363.5 

-46.0 

-43-9 

-45-3 

-45.4 

552.1 

-48.0 

-45.8 

-47.2 

-47.4 

541  .0 

-50.0 

-47.7 

-49.2 

-49.4 

330.2 

-52.0 

-49.6 

-51.2 

-51.4 

319.7 

-54.0 

-5!  .5 

-53.2 

-53.4 

309.5 

-56.0 

-53.4 

-55.2 

-55.4 

299.6 

-58.0 

-55.3 

-57.2 

-57.3 

289.9 

-60.0 

-57.2 

-59.2 

-59.3 

280.5 

-62.0 

-59.1 

-61-2 

-61.3 

271.3 

-64.0 

-61 .0 

-63.1 

-63.3 

262.3 

—66.0 

-62.9 

-65.1 

-65.3 

253.6 

-68.0 

-64.7 

-67.1 

-67.3 

245.1 

-70.0 

-66-6 

-69.1 

-69.3 

236.8 

-72.0 

-68.5 

-71 . 1 

-71.3 

228.7 

-74.0 

-70.3 

-73.1 

-73.2 

220.8 

-76.0 

-72.2 

-75.1 

-75.2 

213.2 

-78.0 

-74.0 

-77.1 

-77.2 

205.7 

-80.0 

-75.8 

-79.1 

-79.2 

198.4 

-82.0 

-77.7 

-81.1 

-81  .2 

191.3 

-84.0 

-79.5 

-83.0 

-83.2 

184.3 

-86.0 

-81 .3 

-85.0 

-85.2 

177.6 

-88.0 

-83.0 

-87.0 

-87.2 

171.0 

-90.0 

-84.8 

-89.0 

-89.2 

164.6 

-92.0 

-86.5 

-91.0 

-91.2 

158.4 

-94.0 

-88.3 

-93.0 

-93.2 

152.4 

-96.0 

-89.9 

-95.0 

-95.1 
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